Spherical symmetry in /(i?)-gravity 



S. Capozziello**,A. Stabilet^ A. Troisi** 
Dipartimento di Scienze Fisiche and INFN, 
Sez. di Napoli, Universita di Napoli "Federico 11", 
Compl. Univ. di Monte S. Angela, 
Edificio G, Via Cinthia, 1-80126 - Napoli, Italy and 
^ Dipartimento di Fisica "E. R. Caianiello", 

Universita' degli Studi di Salerno, 
Via S. Allende, I-84O8I Baroniasi (SA), Italy. 

Spherical symmetry in f{R) gravity is discussed in details considering also the relations with 
the weak field limit. Exact solutions are obtained for constant Ricci curvature scalar and for Ricci 
scalar depending on the radial coordinate. In particular, we discuss how to obtain results which 
' can be consistently compared with General Relativity giving the well known post-Newtonian and 

. post-Minkowskian limits. Furthermore, we implement a perturbation approach to obtain solutions 

' up to the first order starting from spherically symmetric backgrounds. Exact solutions are given for 

' several classes of f{R) theories in both R = constant and R = R{r). 
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1^ ■ I. INTRODUCTION 

^ ' The recent advent of new cosmological precision tests, capable of probing physics at very large redshifts, has changed 
Q^, the modern view of cosmos ruhng out the Cosmological Standard Model based on General Relativity (GR), radiation 
' and baryonic matter. Beside the introduction of dark matter, needed to fit the astrophysical dynamics at galactic 
bJT)! and cluster scales (i.e. to explain clustered structures), a new ingredient is requested in order to explain the observed 
accelerated behavior of the Hubble flow: the dark energy. 
' In particular, the luminosity distance of la Type Supernovae the Large Scale Structure [2|] and the anisotropy 
>^ . of Cosmic Microwave Background [3] suggest that the widely accepted Cosmological Concordance Model (ACDM) is 
^— I ' a spatially flat Universe, dominated by cold dark matter (CDM (~ 0.25 0.3%) which should explain the clustered 
0^ ■ structures) and dark energy (A ~ 0.65 0.7%), in the form of an "effective" cosmological constant, giving rise to the 
00 ' accelerated behavior. 

" Although the cosmological constant [1, H, 0] remains the most relevant candidate to interpret the accelerated behav- 
ior, several proposals have been suggested in the last few years: quintessence models, where the cosmic acceleration 
is generated by means of a scalar field, in a way similar to the early time inflation, acting at large scales and recent 
epochs 0; models based on exotic fluids like the Chaplygin-gas or non-perfect fluids @); phantom fields, based 
on scalar fields with anomalous signature in the kinetic term higher dimensional scenarios (braneworld) 
Actually, all of these models, are based on the peculiar characteristic of introducing new sources into the cosmological 
dynamics, while it would be preferable to develop scenarios consistent with observations without invoking further 
parameters or components non-testable (up to now) at a fundamental level. 
' The resort to modified gravity theories, which extend in some way the GR, allows to pursue this different approach 
(no further unknown sources) giving rise to suitable cosmological models where a late time accelerated expansion 
naturally arises. 

The idea that the Einstein gravity should be extended or corrected at large scales (infrared limit) or at high 
energies (ultraviolet limit) is suggested by several theoretical and observational issues. Quantum field theories in 
curved spacetimes, as well as the low energy limit of string theory, both imply semi - classical effective Lagrangians 
containing higher-order curvature invariants or scalar-tensor terms. In addition, GR has been tested only at solar 
system scales while it shows several shortcomings if checked at higher energies or larger scales. 

Of course modifying the gravitational action asks for several fundamental challenges. These models can exhibit 
instabilities or ghost -like behaviors [Hj, while, on the other side, they should be matched with the low energy 
limit observations and experiments (solar system tests, PPN limit). Despite of all these issues, in the last years, 



o 



- '—I 

X 



* e-mail address: capozziello@na.infn.it 
t e - mail address: stabile@sa.infn.it 
t e-mail address: antrois@gmail.com 



2 



several interesting results have been achieved in the framework of the so called f{R)- gravity at cosmological, galactic 
and solar system scales. 

For example, models based on generic functions of the Ricci scalar R show cosmological solution with late time 
accelerating dynamics ^Tg'], in addition, it has been shown that some of them could agree with CMBR observational 
prescriptions nevertheless this matter is still argument of debate [3]. For a review of the models and their 
cosmological applications see, e.g., 

Moreover, considering /(_R)-gravity in the low energy limit, it is possible to obtain corrected gravitational potentials 
capable of explaining the flat rotation curves of spiral galaxies without considering huge amounts of dark matter 
[2ll . [2^ m, (ij and, furthermore, this seems the only self-consistent way to reproduce the universal rotation curve of 
spiral galaxies 12511 . On the other hand, several anomalies in solar system experiments could be framed and addressed 
in this picture [26| . 

In the last few years, several authors have dealt with this matter considering the Parameterized Post Newtonian 
(PPN) limit [131 . On the other hand, the investigation of spherically symmetric solutions for such kind of models has 
been developed in several papers [H, [sO, [sU, [321 ■ Such an analysis deserves particular attention since it can allow 
to draw interesting conclusions on the effective modification of the gravitational potential induced by higher order 
gravity at low energies and, in addition, it could shed new light on the PPN limit of such theories. However, in several 
recent papers, it has been shown that several classes of f{R) models fairly evade the Solar System constraints and 
agree with the limits of PPN-parameters, in particular with 7^1 probed by experiments [s^ [s^, [13, [H, [s^ . This is 
not the argument of the present paper, however it is important to stress that the spherical symmetry and the weak 
filed limit have to be carefully considered in order to find out physically viable models. 

In this paper, we are going to analyze, in a general way and without specifying a priori the form of the Lagrangian, 
the relation between the spherical symmetry and the weak field limit of f{R) gravity. Our aim is to develop a 
systematic approach considering the theoretical prescriptions to obtain a correct weak field limit in order to point 
out the analogies and the differences with respect to GR. A fundamental issue is to recover the asymptotically flat 
solution in absence of gravity and the well-known results related to the specific case f{R) = R, i.e. GR. Only in this 
situation a correct comparison between GR and any extended gravity theory is possible from an experimental and 
a theoretical viewpoint. For example, as already shown in {m] and discussed in [401, the Birkhoff theorem is not a 
general result for fourth-order models also if it holds for several interesting classes of these theories as discussed, for 
example, in [H, [HI, [H . 

The layout of the paper is the following. Sec. II is devoted to some general remarks on spherical symmetry in f{R) 
gravity. In particular, we derive the field equations in such a symmetry. In Sec. Ill, the expression of the Ricci scalar 
and the general form of metric components are derived in spherical symmetry discussing how recovering the correct 
Minkowski flat limit. Sec. IV is devoted to the discussion of spherically symmetric background solutions with constant 
scalar curvature considering, in particular Schwarzschild-like and Schwarzschild-de Sitter-like solutions with constant 
curvature. Some remarkable /(i?) models are worked out as example. In Sec. V, we discuss the cases in which the 
spherical symmetry is present also for the Ricci scalar depending on the radial coordinate r. This is an interesting 
situation, not present in GR. In fact, as it is well known, in the Einstein theory, the Birkhoff theorem states that 
a spherically symmetric solution is always stationary and static [44J and the Ricci scalar is constant. In f{R) the 
situation is more general and then the Ricci scalar, in principle, can evolve with radial and time coordinates. Sec. VI 
is devoted to the study of a perturbation approach starting from a spherically symmetric background considering the 
general case in which the Ricci scalar is R — R{r). The motivation is due to the fact that, in GR, the Schwarzschild 
solution and the weak held limit coincide under suitable conditions. This could be a good test bed to construct 
a well-posed weak fleld limit for any extended theory of gravity. Finally, we work out some /(i?)-models finding 
spherically-symmetric solutions by the perturbation approach. Conclusions are drawn in Sec. VII. 



II. SPHERICAL SYMMETRY 



Let us consider an analytic function f{R) of the Ricci scalar R. The variational principle for this action is: 




(1) 



where X = 



, Cm is the standard matter Lagrangian and g is the determinant of the metric. By varying with 
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respect to the metric, we obtain the field equations ^ 



f'{R)Rf.. - \f{R)9^.. - f'iR);^^ + 9>^Mf'{R) = ^T^u , (2) 



where the trace equation read 



iUfiR) + f'{R)R - 2f{R) = |t , (3) 



-2 Sj^C r, 



and T^i, = —j= ^ is the energy momentum tensor of standard matter, T = T'^ a is the related trace and 



f'{R) = ^^^w'- Eqs-lH]) can be rewritten in an Einstein- like form recasting the higher than second order contributions 



dR , 

under the form of an effective stress - energy tensor of geometrical origin 



G,. = R,. \r9,. = rir™^ + T^ir^ • (4) 
The effective curvature stress - energy tensor is [l^ : 

r^r"^ = {9,. [f{R) - Rfm + nRY'" {9,p9ua - 9pag,.)} ■ (5) 

The matter term enters Eqs.Q through the modified stress - energy tensor: 

2/'(i?) ' ^ ' 

where there is a non-minimal coupling between matter and geometric degrees of freedom. For the purposes of this 
paper, we will refer to the field equations in the form ([2]) which reveal to be more useful for the following considerations. 
The most general spherically symmetric metric can be written as follows : 

ds^ = mi{t\r')dt''^ +rn2{t' ,r')dr''^ +m3,{t' ,r')dt'dr' +mi{t' ,r')dn, (7) 

where rrii are functions of the radius r' and of the time t' (see also [s^]). In GR, due to the freedom in the frame 
choice, we can consider a coordinate transformation t = Ui{t' , r') , r = U2{t' , r') which maps the metric ([T]) in a new 
one where the off- diagonal term mj,{t' ^r')dt' dr' vanishes and mi{t' ,r') = — t"^, that is^ : 

ds^ = Ait, r)dt^ - B{t, r)dr^ - r^dVt . (8) 

This expression can be considered without loss of generality as the most general definition of a spherically symmetric 
metric compatible with a pseudo - Riemannian manifold without torsion. Actually, by inserting this metric into the 
field equations ([2]) and ([3]), one obtains: 

- \f{R)g^u + n^. = jT^. , (9) 
g-'^H^r = f'iR)R - 2/(i?) ^T, (10) 



^ It is possible to take into account also the Palatini approach in which the metric g and the connection F are considered independent 
variables (see for example |45|I '). Here we will consider the Levi-Civita connection and will use the metric approach. See Il9l |46|| for a 
detailed comparison between the two pictures. 

^ This condition allows to obtain the standard definition of the circumference with radius r. 
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where the two quantities 'H^v and Ti read : 



T^tiu = -f"{R){ R.iiu ~ r* i?^f - ^''uvR,r - gfj.u 



/'"(i?) 



(11) 



n = g^^n^r = 3/"(i?) 



-3/"'(i?) 



(12) 



where V^^ are the standard Christoffel's symbols related to the metric g^i/. These equations show the interesting 
feature that the derivatives with respect to R of f{R) are very well distinct with respect to the time and spatial 
derivatives of R. This peculiarity will allow to better understand the dynamical behavior of solutions, as we shall see 
below. 



III. THE RICCI CURVATURE SCALAR IN SPHERICAL SYMMETRY 



As standard, the Ricci scalar can be written as 

R = 9'"'R^.. - g''-' 

where F"^ is the Christoffel symbols defined as 

1 



■pa -pQ I p/3 pa pa p/3 

/ii/,a fia.iy a/3 ^i' ^ol 



Imposing the spherical symmetry ([8]), the Ricci scalar in terms of the gravitational potentials reads: 



(13) 



(14) 



B 



Rit, r) 



AB - A'2 



r{B^ - A'B') + 2B{2A' + rA" - rB) 



4^2 



B^~B + rB' 



2r^A^B^ 



(15) 



where, for sake of brevity, we have discarded the explicit dependence in A{t, r) and B{t, r) and the prime indicates the 
derivative with respect to r while the dot with respect to t. If the metric ([5]) is time- independent, i.e. A{t,r) = a(r), 
B{t,r) = b{r), Ea.lfTS]) assumes the simpler form: 



a{r) 



R{r) 



2b{r) 2a' {r) + ra"{r) - ra'{r)b'{r) 



6(r)a'(r)-^H - 4a(r)^ b{ry - b{r) + rb'{r) 



2r^a{r)^b{ry 



(16) 



where the radial dependence of the gravitational potentials is now explicitly shown. This expression can be seen as a 
constraint for the functions a{r) and b(r) once a specific form of Ricci scalar is given. In particular, it reduces to a 
Bernoulli equation of index two, that is 



b'{r) + h{r)b{r) +l{r)b{rf =0, 



with respect to the metric potential b{r) 



^ l^ r^a'jrf - iajr)^ -2ra{r)[2a{ry +ra{r)"] \^^^^ , f Mr) 



ra{r)['ia{r) +ra'{r)] 



'■R{r) 



4a(r) + ra'{r) 



b{rf = 0. 



(17) 
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A general solution of (fT7|) is: 



&(r) 



exp[— J dr h(r)] 



K + J drl{r) exp[-/ drh{r)] ' 



(18) 



where K is an integration constant while h{r) and l{r) are the two functions which, respectively, define the coefficients 

of the quadratic and the linear term with respect to h{r), as in the standard definition of the Bernoulli equation [47| . 

Looking at the equation, we can notice that it is possible to have l{r) — which implies to find out solutions with a 

2 

Ricci scalar scaling as — ^ term of the radial coordinate. On the other side, it is not possible to have h{r) = 

since otherwise we will get imaginary solutions. A particular consideration deserves the limit r oo. In order to 
achieve a gravitational potential b{r) with the correct Minkowski limit, both h{r) and l{r) have to go to zero provided 
that the quantity r'^R(r) turns out to be constant: this result implies b'{r) — 0, and, finally, also the metric potential 
b{r) has a correct Minkowski limit. 

In general, if we ask for the asymptotic flatness of the metric as a feature of the theory, the Ricci scalar has to 
evolve at infinity as r^" with n ^ 2. Formally it has to be: 



lim r^R{r) 



(19) 



with n € N. Any other behavior of the Ricci scalar could affect the requirement of the correct asymptotic flatness. 
This result can be easily deduced from Eg. p?)) . In fact, let us consider the simplest spherically symmetric case: 



ds' 



a{r)dt^ - 



dr'^ 
air) 



r^dn. 



(20) 



The Bernoulli Eg. p?)) is easy to integrate and the most general metric potential a{r), compatible with the Ricci scalar 
constraint , is : 



a{r) = 1 + — + ^ + ^ 



r^R{r)dr 



dr 



(21) 



where ki and ^2 are integration constants. Actually one gets the standard result a(r) — 1 (Minkowski) for r — > oo 
only if the condition (fT9|) is satisfied, otherwise we get a diverging gravitational potential. 



IV. SOLUTIONS WITH CONSTANT CURVATURE SCALAR 

The case of constant curvature is equivalent to GR with a cosmological constant and the solution is time independent. 
This result is well known (see, for example, ^A^) but we report, for the sake of completeness, some considerations 
related with it in order to deal with more general case where a radial dependence is supposed. Let us assume a scalar 
curvature constant (i? — Rq). The field Eqs.® and (|10p . being = 0, reduce to: 



foRp.l' r,fo9fJ.l' 



(22) 



X 

foRo - 2/o = yT. 

where /(i?o) — fo, f'{Ro) = fo- Such equations can be arranged as: 



(23) 



X 



Ro = q^T - 4A 



(24) 
(25) 
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where A = ^ttft ^-^d Q ^ ^ fo- standard, the stress-energy tensor of perfect-fluid matter is 
2/0 

Tfj.1, = {p + p)uf_cUu - pgt_Lu ; (26) 

where p is the energy density, p is the pressure and — dx^ / ds is the 4-velocity. A general solution of the above set 
of equations is achieved for p = —p and reads : 



, ^ qXp-2X^,\ dr _ 

\ r 6 / g.yp-2A ^2 ^ ' 

This result means that any /(i?)- theory in the case of constant curvature scalar (R — Rq) exhibits solutions with 
cosmological constant - like behavior (de Sitter) . This is one of the reason wh y th e dark energy issue can be addressed 
using these theory [l^, [l^ . This fact is well known using the FRW metric |48[ . As another remark we have to say 
that this solution cannot describe stellar structures since this kind of equation of state does not work for stars so it 
could be interesting only in cosmological contexts. 
If f{R) is analytic, it is: 

f{R) = A + «'oi? + *(i?). (28) 

where 5*0 is a coupling constant, A plays the role of the cosmological constant and ^'(i?) is a generic analytic function 
of R satisfying the condition 



lim R-^^{R) = (29) 

R — ^0 

where ^E*! is a constant. If we neglect the cosmological constant A and ^Pq is set to zero, we obtain a new class of 
theories which, in the limit R ^ 0, does not reproduce GR (from Ea. (P5)) . we have lim/j^o f{R) ~ R^)- In such a case 
analyzing the whole set of Eqs. f^^ and one can observe that both zero and constant 7^ curvature solutions are 
possible. In particular, if i? = i?o = field equations are solved for every form of gravitational potentials entering 
the spherically symmetric background, provided that the Bernoulli equation (jl7p . relating these functions, is fulfilled 
for the particular case R{r) — 0. The solutions are thus defined by the relation 



. , cxp[— r dr h(r)] 

J r[a(r)+ra' (r)] 

being B{t,r) — b{r) from Eq.®. In Table IH we give some examples of /(i?)-theories admitting solutions with 
constantf^ or null scalar curvature. Each model admits Schwarzschild, Schwarzschild- de Sitter, and the class of 
solutions given by ([50]) . 



V. SOLUTIONS WITH CURVATURE SCALAR FUNCTION OF r 



Up to now we have discussed the behavior of f{R) gravity seeking for spherically symmetric solutions with constant 
scalar curvature. This situation is well known in GR and give rise to the Schwarzschild solution (i? = 0) and the 
Schwarzschild -de Sitter solution {R. = Rq ^ 0). The problem can be generalized in f{R) gravity considering the Ricci 
scalar as an arbitrary function of the radial coordinate r. 

This approach is interesting since, in general. Higher Order Gravity theories are supposed to admit such kind of 
solutions and several examples have been found in literature [13, HH, [SO, H^l . Here we want to face the problem from 
general point of view. 

If we choose the Ricci scalar i? as a generic function of the radial coordinate (i? = Ri^)), it is possible to show 
that also in this case the solution of the field Eqs.Q and ^TU\i is time independent (if T^i^ = 0)- In other words, the 
Birkhoff theorem has to hold. The crucial point of the approach is to study the off- diagonal {t, r} - component of Q 
as well as in the case of GR. This equation, for a generic /(i?) reads : 



±(^r'f'{R)^Bit,r) = 0, 



(31) 
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/(i?)- theory: 
R 



Ru.v — 



Field equations: 



i?a^ = 0, with R = 



with R = 0, Ci / 



-Rn^ + A^p^ = with i? : 
= 



(n-2)e2 

with i? = 0, Ci = 
^Rp.v + Xgp.v = with R — Ro, = 0, n — 2 



-Rfi"^ + -^ffM" = with R = 



(m+2)e2 



R 



+ -^Pm^ = 0, with 7? = i?o SO that ^7?"+' + (2 - 71)^2 i?o + (^^i + 2)^3 = 



R^^ = with R^O 

Rfiv + Ag^i, = with i? = — ^ 



R^lv + \gfj.iy = 0, with R= — ^ 



TABLE I: Emamples of /(_R)-models admitting constant and zero scalar curvature solutions. In the right hand side, the field 
equations are given for each model. The power n, m are natural numbers while are generic real constants. 

and two possibilities are in order. Firstly, we can choose B(t,r) ^ 0. This choice impUes that J'iR) ~ If this is 

the case, the remaining field equation turn out to be not fulfilled and it can be easily recognized that the dynamical 
system encounters a mathematical incompatibility. 

The only possible solution is given by B{t, r) = and then the gravitational potential has to be B{t,r) = b{r). 
Considering also the {6, 0} - equation one can determine that the gravitational potential A{t, r) can be factorized with 
respect to the time, so that we get solutions which can be recast in the stationary spherically symmetric form after a 
suitable coordinate transformation. 

As a matter of fact, even the more general radial dependent case admit time - independent solutions. From the 
trace equation and the {0, 0} - component, we deduce a relation which links A(t, r) — a(r) and B{t, r) — b{r) : 



a{r) = 



b{r)e 



2 r (fi/'-2/)b(r) , 

3 J R'f" 



(32) 



(with /" / 0) and one which relates 6(r) and /(i?) (see also [301 fo'^ ^ similar result) : 



rfirR'f" -Af') + 2f'{rRif' -rR'n-SR'f")- ^ ' 

As above, three further equations has to be satisfied to completely solve the system (respectively the {t, t} and 
{r, r} components of the field equations and the Ricci scalar constraint) while the only unknown functions arc f{R) 
and the Ricci scalar R{r). 

If we now consider a fourth order model of the form f{R) = i?+ $(i?), with $(i?) <C i? we are capable of satisfying 
the whole set of equations up to third order in $. In particular, we can solve the whole set of equations; the relations 
([32|) and (l33|) will give the general solution depending only on the forms of <&(i?) and R — R{r), that is: 
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b{r)e 



2 f [R + (2<!> - m')]b(r) , 
dr 



3 / 



air) = — — ^ — ^ (34) 

3(ri?-<&'0,. 

5(r) = (35) 

This solution is one of the main results of this paper. In fact, once the radial dependence of the scalar curvature is 
obtained, Eqs.([34]) allow to write down the solution of the field equations and the gravitational potential, related to 
the function a(r), can be deduced. Furthermore one can check the physical relevance of such a potential by means of 
astrophysical data, see for example the analysis in [3l|. As a final remark, we have to stress that the solution ([34]) 
has been achieved here in the vacuum case but, in general, R{r) should be obtained from the source T using the 
trace Eq.([3]). This task is simpler in Palatini approach than in the metric approach since, in the former case, the box 
operator is not present in the trace equation. For a detailed discussion, see [28l[29l|. 

VI. PERTURBING THE SPHERICALLY SYMMETRIC SOLUTIONS 

The search for solutions in /(i?)-gravity, in the case of Ricci scalar dependent on the radial coordinate, can be 
faced by means of a perturbation approach. There are several perturbation techniques by which higher order gravity 
can be investigated in the weak field limit. A general approach is starting from analytical f{R) theories assuming 
that the background model slightly deviates from the Einstein GR (this means to consider f{R) = R + <&(-R) where 
$(i?) <C -R as above). Another approach can be developed starting from the background metric considered as the 
0th - order solution. Both these approaches assume the weak field limit of a given higher order gravity theory as a 
correction to the Einstein GR, supposing that zero order approximation should yield the standard lore. 

Both these methods can provide interesting results on the astrophysical scales where spherically symmetric solutions 
characterized by small values of the scalar curvature, can be taken into account. 

In the following, we will consider the first approach assuming that the background metric matches, at zero order, 
the GR solutions. 

In general, searching for solutions by a perturbation technique inGRns to perturb the metric — 9^^^ ^~ ■ 

This 

implies that the field equations ^ and (fTO]) split, up to first order, in two levels. Besides, a perturbation on the 
metric acts also on the Ricci scalar R (see the relation (fT5)) ) and then we can Taylor expand the analytic f{R) about 
the background value of R, i.e.: 



f(^R) = yLlLR^RiO) =y-:^I^(l)". (36) 
n '- n 

However the above condition ^(i?) <C R has to imply the validity of the linear approximation /"(i?(°))//'(i?(°))i?(^^ <C 
1. This is demonstrated by assuming /'(i?) = l+<I>'(i?) and/"(i?) = $"(i?). Immediately we obtain that the condition 
is fulfilled for 

$"(i?("))i?,(i) , , 

For example, given a Lagrangian of the form f{R) = R + —, it means 

R 



"(0) 



j-n{Q) 



2fiR 



(1) 



R(0){R{oy 

while, for /(i?) = i? + ai?2, it is 

2ai?(i' 



< 1 , (38) 



2ai?(0) 



< 1 . (39) 



This means that the validity of the approximation strictly depends on the form of the models and the value of the 
parameters, in the previous case fi and a. For the considerations below, we will assume that it holds. A detailed 
discussion for the Palatini formalism is in [23, Urn ■ 
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The zero order field equations read : 



f'(0) d(0) _ i„(0) .(0) , 0,(0) ^ ^^(0) 



where 



(40) 



-/"'(«) - ^ . (41) 



At first order one lias : 



I'^'Ar^} - + r^i^J^ - i/("\g« + ^1,1^ = f T« (42) 



with 



+ (in V^|>« + In ^/^;>(S)) +5(1)- In ^/^!>(5)] - (^^^^-.i?!^) 

+g^°^^^R%+g(°^''^lnV^g^^R^j)yr^'^R^'^{R(°J^ (43) 



A part the analyticity, no hypothesis has been invoked on the form of f{R). As a matter of fact, f{R) can be 
completely general. At this level, to solve the problem, it is required the zero order solution of Eqs. (|40p which, in 
general, could be a GR solution. This problem can be overcome assuming the same order of perturbation on the 
f{R), that is: 

/ = i? + $(i?), (44) 
where $(-R) is a generic function of the Ricci scalar fulfilling the prescription $ ^ _R. Then we have 

/ - + + , /' = 1 + , /" - , /"' - , (45) 

and Eqs. (|^0)l reduce to the form 

R^ - \R^'^g^$ - jTioj . (46) 



(47) 



On the other hand, Eqs. pzj) reduce to 



where 
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+g^'^>^R^°l + In V^J.^'i?^^ | . (48) 

The new system of field equations is evidently simpler than the starting one and once the zero order solution is 
obtained, the solutions at the first order correction can be easily achieved. In Table IIIl a list of solutions, obtained 
with this perturbation method, is given considering different classes of f{R) models. 

Some remarks on these solutions are in order at this point. In the case of f{R) models which are evidently corrections 
to the Hilbert-Einstein Lagrangian as A + R + eR\n R and R + ei?", with e <C 1, one obtains exact solutions for the 
gravitational potentials a{r) and 6(r) related by a{r) = b(r)~^. The first order expansion is straightforward as in GR. 
If the functions a{r) and 6(r) are not related, for f{R) = A + R + eR\n R, the first order system is directly solved 
without any prescription on the perturbation functions x{r) and y{r). This is not the case for f{R) = R + eR^ since, 
for this model, one can obtain an explicit constraint on the perturbation function implying the possibility to deduce 

2(j){r) 

the form of the gravitational potential 0(r) from air) = 1 H ^ — . In such a case, no corrections are found with 

R 

respect to the standard Newtonian potential. The theories f{R) = i?" and f{R) ~ ^ show similar behaviors. 
The case f{R) — R^ is peculiar and it has to be dealt independently. 
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fiR) - theory: 
spherical potentials: 
solutions: 
first order metric: 

solutions: 

f{R) - theory: 
spherical potentials: 
solutions: 
first order metric: 
solutions: 

f{R) - theory: 
spherical potentials: 

solutions: 

first order metric: 

solutions: 
first order metric: 

solutions: 

fiR) - theory: 
first order metric: 

solutions: 



A + R + eRlnR 
a{r) = b{r)-' = 1 + ^ - ^ + 6x{r) 
x{r) = ^ + ^'^\^-(-lf)-^r- 
a(r) = 1 - ^ + &x{r), b{r) = + Sy{r) 

^ 6 

= C^j,2 _ g^ fj^ , r . 4a(2A^r^-15Ar'' + 18)./(r) + r{36reA[log(-2A)-ll+a(Ar''-6)^./'(r)} 
•'V/ V / ^ ^ ~^ J 36T-<5(Ar-^ —6) 

/ X _ fc2^-6r3eA[ln(-2A)-l] 



a(r) = b(r)-^ = l + ^+5x{r) 
x{r) = !f 
a(r) = l + 5^, b{r) = l+5^ 
x{r) = fci + /c2r, y{r) = fea 

a(r) = b{r)-' = 1 + + Ikf + 5x{r) 

3fc2— fca _i_ k^r'^ I fc4 



n = 2, i?o / and x{r) = ^^^^^ + + ^ / rf,- r^^ f dr- 

with ro satisfying the condition 6ki + 8ro + -Ro^o = 
n'>2, System solved only whit Ro = and no prescriptions on x{r) 



a{r) = l + 6^, bir) = l+6^ 

n = 2 y{r) = - _ ifcl + ^rx'{r) + ki, R{r) = 5Ro 

y{r) = - | / drr^R(r) - ^ + \rx'(r) + fci with R{r) whatever 

a(r) = 1 - ^ + 5x{r), b{r) = ^ + 5j/(r) 



n = 2 2/(r) = 



3r'^-7rgr+4r^ 3(3t-^ -7T-gr-+4T-^) 

\n^2 whatever functions x{r) , y{r) and -R(r') 



Ji/(? + i?) 



(3rg-4r)(r-g-r)^ 



a{r) = l + d^, bir) = l+6^ 



a::(r) 
2/(?') 



4e 

e 

2e V6 (6gl/^ + V6gr) ,^ 
36372 «1 



2V^e \/6 
J372 



^2 + ksr 



{1/2,. 

2e n/8 (VE-i'^^^r) 
5372 



k2 



TABLE II: A list of exact solutions obtained via the perturbation approach for several classes of f(R) theories; ki are integration 
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VII. CONCLUSIONS 



General Relativity has been consistently tested in physical situations implying, essentially, spherical symmetry 
and weak field limit [i^. One of the fundamental and obvious issue that any theory of gravity should satisfy is 
the fact that, in absence of gravitational field or very far from a given distribution of sources, the spacetime has 
to be asymptotically flat (Minkowski). Any alternative or modified gravitational theory (beside the diffeomorphism 
invariance and the general covariance) should address these physical requirements to be consistently compared with 
GR. This is a crucial point which several times is not considered when people is constructing the weak field limit of 
alternative theories of gravity. 

In this paper, we have faced this problem discussing, in the most general way, what the meaning of spherical 
solutions is in f{R) theories of gravity and when the standard results of GR are recovered in the limits r —^ oo and 
f{R) R. Essentially, spherical solutions can be classified, with respect to the Ricci curvature scalar R, a,s R — 0, 
R = Rq ^ 0, and R = R{r), where Rq is a constant and R{r) is a function of the radial coordinate r. In these cases, 
the Birkhoff theorem holds; this means that stationary solutions are also static. However, as shown in the companion 
paper [s^ . this theorem does not hold, in general, for every f{R) theory since time-dependent evolution can emerge 
depending on the order of perturbations. 

In order to achieve exact spherical solutions, a crucial role is played by the relations existing between the metric 
potentials and between them and the Ricci curvature scalar. In particular, the relations between the metric potentials 
and the Ricci scalar can be used as a constraint: this gives a Bernoulli equation. Solving it, in principle, spherically 
symmetric solutions can be obtained for any analytic f{R) function, both for constant curvature scalar and for 
curvature scalar depending on r. 

Such spherically symmetric solutions can be used as background to test how generic f{R) theories of gravity deviate 
from GR. Particularly interesting are those theories that imply f{R) R in the weak field limit. In such cases, the 
experimental comparison is straightforward and also experimental results, evading GR constraints, can be framed in 
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Finally, we have constructed a perturbation approach in which we search for spherically symmetric solutions at 
the Oth-order and then we search for solutions at the first order. The scheme is iterative and could be, in principle, 
extended to any order in perturbations. The crucial request is to take into account f{R) theories which are Taylor 
expandable about some value R — Rq of the curvature scalar. 

A important remark is in order at this point. Considering interior and exterior solutions, the junction conditions are 
related to the integration constants of the problem and strictly depend on the source (e.g. the form of T). We have 
not considered this aspect here since we have, essentially, searched for vacuum solutions. However, such a problem 
has to be carefully faced in order to deal with physically consistent solutions. For example, the Schwarzschild solution 
R = 0, which is one of the exterior solutions which we have considered, always satisfies the junction conditions with 
physically interesting interior metric. This is not the case for several spherically symmetric solutions which could give 
rise to unphysical junction conditions and not be in agreement with Newton's law of gravitation, also asymptotically. 
In these cases, such solutions have to be discarded. A detailed analysis in this sense is in the paper for spherically 
symmetric solutions obtained in the Palatini formalism and in [Zsj for the metric approach. 
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